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Abstract 

An angular analysis of the B° —> K*°e + e~ decay is performed using a data sample, 
corresponding to an integrated luminosity of 3.0 fb _1 , collected by the LHCb 
experiment in pp collisions at centre-of-mass energies of 7 and 8 TeV during 2011 
and 2012. For the first time several observables are measured in the dielectron mass 
squared ( q 2 ) interval between 0.002 and 1.120 GeV 2 /c 4 . The angular observables Fl 
and A^ c which are related to the K*° polarisation and to the lepton forward-backward 
asymmetry, are measured to be El = 0.16 ±0.06 ±0.03 and A^ e = 0.10 ±0.18 ±0.05, 
where the first uncertainty is statistical and the second systematic. The angular 
observables Aij 1 and which are sensitive to the photon polarisation in this q 2 

range, are found to be = -0.23 ± 0.23 ± 0.05 and A\ n = 0.14 ± 0.22 ± 0.05. The 
results are consistent with Standard Model predictions. 
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1 Introduction 


The B° —> K*°e + e~ decay is a flavour changing neutral current process that is mediated 
by electroweak box and loop diagrams in the Standard Model (SM). Charge conjugation is 
implied throughout this paper unless stated otherwise and the K*° represents the Ji*°(892), 
reconstructed as K*° —>• K + tt~. The angular distribution of the K + n~e + e~ system is 
particularly sensitive to contributions from non-SM physics (NP). The leading SM diagrams 
are shown in Fig. [lj the relative contribution of each of the diagrams varies with the 
dilepton invariant mass. In the region where the dilepton invariant mass squared (g 2 ) is less 
than 6 GeV 2 /c 4 , some theoretical uncertainties from long distance contributions are greatly 
reduced, thereby allowing more control over the SM prediction and increasing sensitivity 
to any NP effect |1,2|. Furthermore, the contribution from a virtual photon coupling 
to the lepton pair dominates in the very low g 2 region, allowing measurement of the 
helicity of the photon in b—> s'y transitions [3.4 . In the SM, this photon is predominantly 
left-handed, with a small right-handed component arising from the mass of the s quark 
and long distance effects. In contrast, in many extension of the SM, NP may manifest as 
a large right handed current, see for example Refs. [5j[8]. 

The q 2 region below lGeV 2 /c 4 has previously been studied through the analysis of 
the -B—)• K*£ + e~ (£ = e,n) &ui- Experimentally, an analysis with muons rather than 
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Figure 1: Dominant Standard Model Feynman graphs for the electroweak loop and box diagrams 
involved in the B° —> K*°e + e~ decay. 
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electrons in the final state produces a much higher yield at LHCb. This is primarily 
due to the distinctive signature that muons provide, which is efficiently exploited in 
the online selection, together with the better mass and energy resolutions and higher 
reconstruction efficiency of dimuon decays. However, as outlined in Ref. 12 , dielectron 


decays at low q 2 provide greater sensitivity to the photon polarisation and therefore to the 
C 7 and C 7 Wilson coefficients, which are associated with the left-handed and right-handed 
electromagnetic operators, respectively |3|. Due to the muon mass, the virtual photon 
contribution in dimuon decays is suppressed relative to dielectron decays. Additionally, the 
formalism for the B° —)■ K*°e + e~ decay is greatly simplified as the electron mass can be 
neglected. Indeed, the decay with electrons allows for an angular analysis down to a q 2 of 
0.0004 GeV 2 /c 4 . However, above a q 2 of 1 GeV 2 /c 4 , the muon mass terms become negligible 
and the electron and muon modes have essentially the same functional dependence on the 
Wilson coefficients (within the lepton flavour universality assumption). 

This work is based on a previous analysis performed by the LHCb collaboration to 
measure the B° —>• K*°e + e~ branching fraction with an integrated luminosity of 1.0 fb _1 13 , 
with the selection re-optimised for the angular analysis. 

The partial decay width of the B° —> K*°e + e~ decay can be described in terms of q 2 
and three angles, 0£, 9 K and 0. The angle 9g is defined as the angle between the direction 
of the e + (e _ ) and the direction opposite to that of the B° ( B° ) meson in the dielectron 
rest frame. The angle 9k is defined as the angle between the direction of the kaon and the 
direction opposite to that of the B° ( B° ) meson in the K*° (K*°) rest frame. The angle 
0 is the angle between the plane containing the e + and e~ and the plane containing the 
kaon and pion from the K*° (K*°) in the B° ( B° ) rest frame. The basis is designed such 
that the angular definition for the B° decay is a CP transformation of that for the B° 
decay. These definitions are identical to those used for the R° —s- K* 0 /j, + /j,~ analysis 11 
As in Ref. 
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the angle 0 is transformed such that 0 = 0 + n if 0 < 0, to compensate 
for the limited signal yield. This transformation cancels out the terms that have a sin0 
or cos 0 dependence and simplifies the angular expression without any loss of sensitivity 
to the remaining observables. In the limit of massless leptons and neglecting the K + n~ 
S-wave contribution, which is expected to be negligible^ at low q 2 with the current sample 


1 Using Refs |T 
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it can be shown that the ratio of the S-wave fraction to the fraction of longitudinal 
polarisation of the K*° is constant as function of q 2 in the 0-6 GeV 2 /c 4 range. 
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14 , the B° —> K*°e + e angular distribution reads as 

1 d 4 (r + f) 
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d(r + T)/dg 2 dg 2 dcos 9g dcos 9k d 0 16?r 


-(1 - F L ) sin 2 9 k + F l cos 2 9 k + 

^ (1 — F l ) sin 2 9 k — F L cos 2 9 K ^j cos 29 e + 

^(1 — F l )A^ sin 2 0 K sin 2 9 1 cos 20 + 

(1 — F l ) A^ e sin 2 9k cos 9g + 


-(1 — Fl)A^ sin 2 9k sin 2 9g sin 20 


( 1 ) 

The four angular observables Fl, A^\ and A are related to the transversity 
amplitudes through |2 


F _ l^ol 2 _ 

L \Ao\* + \A n \* + \A ± \* 

A (2) I^-l| 2 ~ l^lll 2 

T l^| 2 +l^||| 2 

4Re _ “FFe{A\\ lA±l + A\\ R A* ±R ) 

T “ l^||l 2 + l^l 2 
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( 2 ) 


where \A 0 \ 2 = \A ol \ 2 + \A or \ 2 , \A ±\ 2 = \A ±L \ 2 + (A^ 2 and |A||| 2 = \A\\ L \ 2 + |zl|| R | 2 . The 
amplitudes Ao, ^4|| and A± correspond to different polarisation states of the K*° in the 
decay. The labels L and R refer to the left and right chirality of the dielectron system. 

Given the definition of 0, the observable Aj?' 1 is averaged between B° and B° decays, 
while A Ip corresponds to a CP asymmetry 1151. The observable Fl is the longitudinal 
polarisation of the K*° and is expected to be small at low q 2 , since the virtual photon 
is then quasi-real and therefore transversely polarised. The observable is related to 
the forward-backward asymmetry A FB by = |A fb /(1 — F L ) [2|. The observables 
and Alp, in the limit q 2 —> 0, can be expressed as simple functions of the C 7 and C 7 
coefficients |2 


4V 0) 


21Ze(C 7 C 7 *) 


and A^(q 2 ->• 0) 


21rn(C 7 C 7 *) 

ic 7 | 2 + ic;i 2 


(3) 


These measurements therefore provide information on photon polarisation amplitudes, 
similar to that obtained by the CP asymmetry measured through time-dependent analyses 
in B°— » K*°(—> Jig7r°)7 decays 

This paper presents measurements of Fl, A^\ A If 1 and Alp 6 of the B° —» K*°e + e~ 
decay in the bin corresponding to a reconstructed q 2 from 0.0004 to 1 GeV 2 /c 4 . 
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2 The LHCb detector and data set 


The study reported here is based on pp collision data, corresponding to an integrated 
luminosity of 3.0 fb -1 , collected at the Large Hadron Collider (LHC) with the LHCb 
detector 18,19 at centre-of-mass energies of 7 and 8 TeV during 2011 and 2012. The 


LHCb detector is a single-arm forward spectrometer covering the pseudorapidity range 
2 < t] < 5, designed for the study of particles containing b or c quarks. The detector includes 
a high-precision tracking system consisting of a silicon-strip vertex detector surrounding 
the pp interaction region [20 , a large-area silicon-strip detector located upstream of a 


dipole magnet with a bending power of about 4 Tin, and three stations of silicon-strip 
detectors and straw drift tubes |2l] placed downstream of the magnet. The tracking system 
provides a measurement of momentum, p, with a relative uncertainty that varies from 0.5% 
at low momentum to 1.0% at 200 GeV/c. The minimum distance of a track to a primary 
vertex, the impact parameter (IP), is measured with a resolution of (15 + 29/px) M-nr, where 
Pt is the component of the momentum transverse to the beam, in GeV/c. Different types 
of charged hadrons are distinguished using information from two ring-imaging Cherenkov 
detectors 22 . Photons, electrons and hadrons are identified by a calorimeter system 


consisting of scintillating-pad and preshower detectors, an electromagnetic calorimeter 
(ECAL) and a hadronic calorimeter. Muons are identified by a system composed of 
alternating layers of iron and multiwire proportional chambers 23 . 

The trigger [24 consists of a hardware stage, based on information from the calorimeter 


and muon systems, followed by a software stage, which applies a full event reconstruction. 
For signal candidates to be considered in this analysis, all tracks from the B° —> K*°e + e~ 
decay must have hits in the vertex detector and at least one of the tracks from the 
B° —> K*°e + e~ decay must meet the requirements of the hardware electron or hadron 
triggers, or the hardware trigger must be fulfilled independently of any of the decay 
products of the signal B° candidate (usually triggering on the other b hadron in the 
event). The hardware electron trigger requires the presence of an ECAL cluster with 
a minimum transverse energy between 2.5 GeV and 2.96 GeV depending on the data 
taking period. The hardware hadron trigger requires the presence of a cluster in the 
hadron calorimeter with a transverse energy greater than 3.5 GeV. The software trigger 
requires a two-, three- or four-track secondary vertex with a significant displacement from 
the primary pp interaction vertices (PVs). At least one charged particle must have a 
transverse momentum px > 1-7 GeV/c and be inconsistent with originating from the PV. 
A multivariate algorithm 25 is used for the identification of secondary vertices consistent 


with the decay of a b hadron. 

Samples of simulated B° —> K*°e + e~ events are used to determine the efficiency to 
trigger, reconstruct and select signal events. In addition, specific samples of simulated 
events are utilised to estimate the contribution from exclusive backgrounds and to model 


28 


their mass and angular distributions. The pp collisions are generated using PYTHIA 26 
with a specific LHCb configuration 
by EvtGen 


27]. Decays of hadronic particles are described 
in which final-state radiation is generated using PHOTOS 


29 . The 


interaction of the generated particles with the detector, and its response, are implemented 
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using the Geant4 toolkit 1301 as described in Ref. 31 . The simulated samples are 


corrected for known differences between data and simulation in particle identification 
detector occupancy and hardware trigger efficiency. 
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3 Selection of signal candidates 

Bremsstrahlung radiation, if not accounted for, would worsen the B° mass resolution. If 
the radiation occurs downstream of the dipole magnet, the momentum of the electron is 
correctly measured and the photon energy is deposited in the same calorimeter cell as 
the electron. If photons are emitted upstream of the magnet, the electron momentum is 
evaluated after photon emission, and the measured B° mass is shifted. In general, these 
bremsstrahlung photons deposit their energy in different calorimeter cells than those hit by 
the electron. In both cases, the ratio of the energy detected in the ECAL to the momentum 
measured by the tracking system, an important variable in identifying electrons, remains 
unbiased. To improve the momentum reconstruction, a dedicated bremsstrahlung recovery 
is used. Contributions from photon candidates, neutral clusters with transverse energy 
greater than 75 MeV, found within a region of the ECAL defined by the extrapolation of 
the electron track upstream of the magnet, are added to the measured electron momentum. 

Oppositely charged electron pairs formed from tracks with px exceeding 350MeV/c and 
with a good-quality vertex are used to form signal candidates. If the same bremsstrahlung 
photon is associated with both the e + and the e~, its energy is added randomly to one of 
the tracks. The reconstructed e + e _ invariant mass is required to be in the range 20-1000 
MeV/c 2 (0.0004 < q 2 < 1 GeV 2 /c 4 ). The choice of the lower bound is a compromise between 
the gain in sensitivity to the photon polarisation from measuring as low as possible in 
q 2 and a degradation of the resolution in <j) as q 2 decreases, due to multiple scattering, 
as shown in Fig. |2j The lower bound requirement at 20 MeV/c 2 on the e + e^ invariant 
mass also serves to reduce the background from B° —>■ K*° 7 decays followed by a photon 
conversion in the material, noted below as B° — » K*° 7 e + e -. 

Candidate K*° mesons are reconstructed in the K*° —s- K + tt~ mode where the px of 
the K + (7T“) meson is required to be larger than 400 (300) MeV/c and charged pions and 
kaons are identified using information from the RICH detectors. 

Candidate K*° mesons and e + e _ pairs are required to have a common good-quality 
vertex to form B° candidates. When more than one PV is reconstructed, the one giving 
the smallest IP y 2 for the B° candidate is chosen. The reconstructed decay vertex of the 
B° candidate is required to be significantly separated from the PV and the candidate 
momentum direction to be consistent with its direction of flight from the PV. The B° mass 
resolution, the angular acceptance and the rates of physics and combinatorial backgrounds 
depend on how the event was triggered. The data sample is therefore divided into three 
mutually exclusive categories: events for which one of the electrons from the B° decay 
satisfies the hardware electron trigger, events for which one of the hadrons from the B° 
decay satisfies the hardware hadron trigger and events triggered by activity in the event 
not due to any of the signal decay particles. 
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Figure 2: Resolution on the cj) angle as a function of the e + e _ invariant mass as obtained from 
LHCb simulated events. 


In order to maximise the signal efficiency while reducing the high level of combina¬ 
torial background, a multivariate classifier based on a boosted decision tree algorithm 
(BDT) 132,33 is used. The signal training sample is composed of simulated B° —> K* 0 ' 


e 1 e 


events and the background training sample is taken from the upper invariant mass sideband 
(m(K + 7r~e + e~) > 5600 MeV/c 2 ) of B° —> K*°e + e~ decays reconstructed in half of the data 
sample. Two separate BDTs are used, one each for half of the data sample. They are 
optimised separately and applied to the complementary half of the data in order to avoid 
any potential bias due to the use of the data upper sideband for the background sample. 
The BDT uses information about the event kinematic properties, vertex and track quality, 
IP and px of the tracks, flight distance from the PV as well as information about isola¬ 
tion of the final state particles]^] The selection is optimised to maximise Ns/s/Ns + Nb 
separately for the three trigger categories and the two BDTs through a grid search of 
the set of criteria for the particle identification of the four final state particles and the 
BDT response. The background yield ( Nb ) is extrapolated into the signal range using 
the m{K + Ti~e + e~) distribution outside a ±300 MeV/c 2 window around the known B° 
mass. The expected signal yield (Ns) is obtained using the I ? 0 —y K*°e + e~ simulation and 
the known B° —> K*°e + e~ branching fraction 113], and correcting for data-to-simulation 

- 'i zv*o 


differences in the selection efficiency obtained using the well known B° —> J/^>(e + e~)K 
decay. The efficiency of this requirement on the selected signal is 93% while the background 
is reduced by two orders of magnitude. The expected values for N s /\JN S + N B range 
from 3.9 to 7.5 depending on the trigger category. 


2 The isolation is defined as the number of good two-track vertices that one of the candidate signal 


tracks can make with any other track in the event 34 . 


6 












4 Exclusive and partially reconstructed backgrounds 


Several sources of background are studied using samples of simulated events, corrected to 
reflect the difference in particle identification performances between data and simulation. 

A large non-peaking background comes from the B° —y D~e + u decay, with 
D~ —y e~VK*° which has a combined branching fraction about four orders of magni¬ 
tude larger than that of the signal. In the rare case where both neutrinos have low energies, 
the signal selection is ineffective at rejecting this background which tends to peak towards 
cos 6b ~ 1. In order to avoid any potential bias in the measurement of the A^ 0 parameter, 
a symmetric requirement of | cos 6p\ < 0.8 is applied to suppress this background, resulting 
in a loss of signal of the order of 10%. 

To suppress background from B° s — > 0e + e _ decays, with 0 — » K + K~, where one of 
the kaons is misidentified as a pion, the two-hadron invariant mass computed under the 
K + K~ hypothesis is required to be larger than 1040MeV/c 2 . 

Background from the decay A® —s- pK~e + e~ is suppressed by rejecting events where 
the pion is consistent with being a proton, according to the information from the RICH 
detectors. 


The probability for a decay B° — >■ K * u e + e~ to be misidentified as R u — > K*°e ' e is 
estimated to be 1.1 % using simulated events and this background is therefore neglected. 

Another important source of background comes from the B° —>• K*° 7 decay, where the 
photon converts into an e + e _ pair. In LHCb, approximately 40% of the photons convert 
before reaching the calorimeter, and although only about 10 % are reconstructed as an e + e~ 
pair with hits in the vertex detector, the resulting mass of the B° candidate peaks in the 
signal region. Two very effective criteria for suppressing this background are the minimum 
requirement on the e + e _ invariant mass, m(e + e _ ) > 20MeV/c 2 , and a requirement that 
the uncertainty of the reconstructed z coordinate of the e + e~ pair, a z (e + e~), is less than 
30 mm. These requirements reject more than 99% of simulated B° —> K *°7 events. The 
remaining contamination is estimated by normalising the simulated B° —> K*°^ e + e - to the 
observed yield without the a z (e + e~ ) criterion and requiring the e + e _ invariant mass to be 
lower than 5 MeV/c 2 . The residual contamination from B° —> K *°7 decays is (3.8 ± 1.9)% 
of the signal yield. Part of this background comes from low-mass e + e~ pairs that are 
reconstructed at larger masses due to multiple scattering. The remainder comes from direct 
Bethe-Heitler pair-production at masses larger than 20 MeV/c 2 . To obtain an accurate 
estimate of this component, the Geant4 simulation is reweighted as a function of the true 
e + e _ mass to match the distribution of Ref. 35 since Geant4 does not model correctly 
the high-mass e + e~ pair production. 

Another possible source of contamination is the decay B° —> K*°V(—> e + e~) where 
V is a p, 00 or 0 meson. Expected rates for these backgrounds have been evaluated in 
Refs. j4 36 . The effects of direct decays or interference with the signal decay are found to 
be negligible after integrating over the q 2 range. 

Partially reconstructed (PR) backgrounds arising from B° —> K*°e + e~X decays, where 
one or more of the decay products (A) from the B° decay is not reconstructed, are also 
taken into account. These incomplete events are mostly due to decays involving higher 
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K* resonances, hereafter referred to as K**. The decays B° —> K*°r] and B° —> K* 0 7i 0 
are also studied and several cases are considered: the case when the e + e _ pair comes 
from a converted photon in the material, the case when the e + and e~ originate from the 
conversions of the two photons and finally the case of the Dalitz decay of the 7j or the 7r°. 
They contribute about 25% of the PR background in the angular fitting domain. 


5 Fit to the K + ti e + e invariant mass distribution 


In a first step, a mass fit over a wide mass range, from 4300 to 6300 MeV/c 2 , is performed to 
estimate the size of the B° —> K*°e + e~ signal, the combinatorial background and the PR 
background. The fractions of each component are determined from unbinned maximum 
likelihood fits to the mass distributions separately for each trigger category. The mass 
distribution of each category is fitted to a sum of probability density functions (PDFs), 
modelling the different components. Following the strategy of Ref. 37 , the signal PDF 


depends on the number of neutral clusters that are added to the dielectron candidate 
to correct for the effects of bremsstrahlung. The signal is described by the sum of a 
Crystal Ball function 1381 (CB) and a wide Gaussian function accounting for the cases 
where background photons have been associated; the CB function accounts for over 90% 
of the total signal PDF. The shape of the combinatorial background is parameterised 
by an exponential function. Finally, the shape of the PR background is described by 


non-parametric PDFs 39 determined from fully simulated events passing the selection. 


The signal shape parameters are fixed to the values obtained from fits to simulation 
but the widths and mean values are corrected for data simulation differences using 
B°J/ip (e + e~)K *° as a control channel. Since the photon pole contribution dominates 
in the low-g 2 region, the PR background is expected to be similar for B° —y K*°e + e~ and 
B° —> K*° 7. The large branching fraction of the decay B° —y K *°7 allows the fractions 
of PR background relative to the signal yield to be determined from the data. These 
fractions are extracted from a fit to a larger sample of events obtained by removing the 
requirements on the lower bound of the e + e~ invariant mass and on cr^(e + e _ ) and therefore 
dominated by B° —> K*° 7 e + e - events. The invariant mass distribution, together with the 
PDFs resulting from this fit, is shown in Fig. |3](a) for the three trigger categories grouped 
together. The corresponding distribution for the R° —> K*°e + e~ fit is shown in Fig. [3](b). 
There are 150 ± 17 B° —> K*°e + e~ signal events, 106 ± 16 PR background events and 
681 ± 32 combinatorial background events in the 4300—6300 MeV/c 2 window. 

In this wide mass window, the sample is dominated by combinatorial background, 
whose angular shape is difficult to model. Furthermore the angular distributions depend 
on the kinematic properties of the background and may thus vary as functions of mass. 
Hence, the angular fit is performed in a narrower mass window from 4800 MeV/c 2 to 
5400 MeV/c 2 . In this restricted window there are 124 B° —y K*°e + e~ signal events, 38 PR 
and 83 combinatorial background events, corresponding to a signal-to-background ratio 
of the order of one. About half of these events belong to the electron hardware trigger 
category and the rest are equally distributed between the other two categories. 
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Figure 3: Invariant mass distribution for (a) the B° —> K*° 7 e + e - and (b) the B° —> K*°e + e~ 
decay modes and the three trigger categories grouped together. The dashed line is the signal 
PDF, the light grey area corresponds to the combinatorial background and the dark grey area 
is the PR background. The solid line is the total PDF. The two vertical dotted lines on the 
B° —> I\*°e + e~ plot indicate the signal window that is used in the angular fit. 


6 Angular acceptance and angular modelling of the 
backgrounds 

6.1 Angular acceptance 

The angular acceptance is factorised as £(cos 9i, cos 9 k, (ft) = £(cos9e)£(cos9K)£((ft) as 
supported by simulation studies. The three corresponding one-dimensional angular dis¬ 
tributions for the B° —s- A'*°e + e~ decay are distorted by the geometrical acceptance of 
the detector, the trigger, the event reconstruction and the selection. Furthermore, their 
precise shapes depend upon the various trigger categories, each being enriched in events 
with different kinematic properties. For the (ft angle, a uniform acceptance is expected. 
However, there are distortions in both the cos 9a and cos 9k distributions, mainly arising 
from requirements on the transverse momenta of the particles. The cos 9k acceptance 
is asymmetric due to the momentum imbalance between the kaon and the pion from 
the K*° decay in the laboratory frame due to their different masses. The cos 9k and 
cos 9e acceptance distributions are modelled on simulated B° —s- K*°e + e~ events with 
Legendre polynomials of fourth order. The functions chosen to model the cos 6^ acceptance 
are assumed to be symmetric and modified by a linear term to estimate the systematic 
uncertainty on the parameter. For the (ft acceptance, no significant deviation from 
uniformity is observed. To estimate the systematic uncertainty, modulations in cos 20 
or sin 20 are allowed. Such modulations are the most harmful ones since they may be 
confused with physics processes yielding non-zero values of or A^ 1 . 
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6.2 Angular modelling of the backgrounds 

In the mass window 4800 < 7n(K + 7i~e + e~) < 5400MeV/c 2 used in the angular analysis, 
about one third of the events are combinatorial background. The angular distribution of 
these events is described by the product of three independent distributions for cos 9g, cos 6 % 
and (j). This background largely dominates at low m{K + Ti~e + e~)\ between 4300MeV/c 2 
and 4800MeV/c 2 , about 90% of the events are combinatorial background according to the 
mass fit shown in Fig. [3j However, the angular distributions of the background depend 
upon m(K + TT~e + e~ ) and the information from the lower mass window cannot be used 
directly for modelling the signal region. The effect of this correlation is extracted from a 
sample of data events selected with a looser BDT requirement but excluding the region 
of the BDT response corresponding to the signal. With this selection the sample is 
dominated by background in the whole mass range. The cos 6 % background distributions 
are modelled as first order polynomials. The cos Qn background distributions are modelled 
with polynomial functions with third and fourth order terms. The (f) distributions are 
compatible with being uniform. This method assumes that there is no strong correlation 
between the BDT response and m{K + ti~ e + e~). This assumption is tested by subdividing 
the sample of events with looser BDT response and comparing the differences between 
the angular shapes predicted by this procedure and those observed. These differences are 
smaller than the statistical uncertainties of the parameters used to describe the angular 
shapes. The statistical uncertainties are thus used to assess the size of the systematic 
uncertainties due to the combinatorial background modelling. 

The PR background accounts for about 15% of the events in the angular fit mass 
window. These events cannot be treated in the same way as the combinatorial ones. Since 
only one or two particles are not reconstructed, the observed angular distributions retain 
some of the features induced by the dynamics of the decay. Hence, they are modelled 
using the same functional shapes as the signal, but with independent physics parameters, 
_Fl,pr, Rx pr; ^t'pr and ^T e pR- The loss of one or more final-state hadrons leads to a 
smaller apparent polarisation of the K*°. While on B° —> K*° 7 simulated events the iy. 
parameter is found to be zero, it reaches 17% for simulated B—t 7/1 **(—)■ KnX) events. 
Since in the SM one expects an F l value of the order of 15 to 20%, Fl,pr is assumed to be 
equal to 1/3, which is equivalent to no polarisation. This parameter is varied between 17% 
and 50% to assess the size of the systematic uncertainty associated with this hypothesis. 
Similarly, the loss of information due to the unreconstructed particles leads to a damping 
of the transverse asymmetries of the PR background, Ht,pr, ^t'pr and Hj/pr) compared 
to those of the signal. The signal transverse asymmetries are expected to be small in 
the SM, therefore their values are set to zero to describe the angular shape of the PR 
background. For pr an d ^t/pr the validity of this assumption is tested by comparing 
angular fits to B —> J/ip K** (—>■ KttX) and B° —> J/t/jK* 0 simulated events, which confirms 
a damping factor compatible with zero. The systematic uncertainty associated with this 
assumption is estimated by varying pr and ^tpr up to half of the fitted signal values 

of and H!,? 1 , i.e. assuming a damping factor of 0.5. For the R^pr parameter, however, 
one cannot estimate a damping factor with the same method since in the B —» J/i/jK* 0 
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decay the value of A3jP is zero. The systematic uncertainty is evaluated by allowing the 
/UjbpR parameter to be as high as the value obtained from the B° —> K*°e + e~ angular 
fit. 

7 Measurement of the angular observables 

7.1 Fit results 

To measure the four angular observables, iy, A^\ and A^ e , an unbinned maximum 
likelihood fit is performed on the m(K + 7i~e + e~), cos 6 b, cos 6k and 0 distributions in 
the signal window defined in Sec. [5] The inclusion of m(K + TT~e + e~) in the fit strongly 
improves its statistical power since the level of background varies significantly within the 
signal mass window. The fit is performed simultaneously on the three trigger categories 
sharing the fit parameters associated with the angular observables. The mass PDFs 
for the three components (signal, PR background and combinatorial background) are 
obtained from the fit described in Sec. [5] The angular PDFs for the signal are obtained 
by multiplying the formula of Eq. [T| by the acceptance described in Sec. [6j Similarly, the 
angular PDFs for the PR background are modelled by using Eq. [I] and the acceptance 
described in Sec. [6 and setting Fl,pr = 0.33 and R^pr = ^t\pr = ^t\pr = 0- Finally, the 
angular PDFs for the combinatorial background are described in Sec. [ 6 ] The combinatorial 
and PR background fractions are constrained to the values calculated from the mass fit 
described in Sec. [5] The fit is validated using a large number of pseudo-experiments that 
include all the components of the fits. Several input values for the angular observables, 
Fl, A^\ and are studied including those associated with NP models, and the fit 
results are in good agreement with the inputs. The fitting procedure is also verified using 
a large sample of fully simulated events; the fitted values of Fl, A^\ A Jp and are in 
excellent agreement with the generated ones. This validates not only the fit but also the 
assumption that the angular acceptance factorises. The distributions of m(K + n~e + e~), 
cos 6 ^, cos 9k and 0, together with the likelihood projections resulting from the fit, are 
shown in Fig. |d]and the fit results are given in Table [lj The fitted values of iR, A^\ 
and are corrected for the (3.8 ± 1.9)% contamination from B° —» K*° 7 e + e - decays, 
assuming that F 0 k*o 7 , ^k*o 7 , ^tk *°7 anc ^ are a ^ e 9 ua l to zero, and are used for 

the computation of the systematic uncertainties related to the angular description of the 
PR background. The fitted values are also corrected for the small fit biases due to the 
limited size of the data sample. 

7.2 Systematic uncertainties 

To evaluate the contributions from the possible sources of systematic uncertainty, pseudo¬ 
experiments with modified parameters are generated and fitted with the PDFs used to 
fit the data. Fit results are then compared with input values to assess the size of the 
uncertainties. 
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Table 1: Fit results for the angular observables i±. and The second column 

corresponds to the uncorrected values directly obtained from the fit while the third column gives 
the final results after the correction for the (3.8 ± 1.9)% of B° —> K* 0r y e + e - contamination and for 
the small fit biases due to the limited size of the data sample. The first uncertainty is statistical 
and the second systematic. 



Uncorrected values 

Corrected values 

F l 

a {2) 

jArji 

4 Im 

ARe 

Jlrji 

0.15 ±0.06 

-0.22 ±0.23 
±0.14 ±0.22 
±0.09 ±0.18 

0.16 ±0.06 ±0.03 

-0.23 ±0.23 ±0.05 
±0.14 ±0.22 ±0.05 
±0.10 ±0.18 ±0.05 
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Figure 4: Distributions of the K + 7T~e + e~ invariant mass, cos 0^, cos 6k and </> variables for the 
B® i^*°e + e _ decay mode and the three trigger categories grouped together. The dashed line 
is the signal PDF, the light grey area corresponds to the combinatorial background, the dark 
grey area is the PR background. The solid line is the total PDF. 


The systematic uncertainties due to the modelling of the angular acceptance are 
estimated by varying the shapes introducing functional dependences that would bias the 
angular observables. 

The uncertainties due to the description of the shape of the combinatorial background 
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Table 2: Summary of the systematic uncertainties. 


Source 


ct(A^) 

a(A!n 

a(A^ e ) 

Acceptance modelling 

0.013 

0.038 

0.035 

0.031 

Combinatorial background 

0.006 

0.030 

0.029 

0.038 

PR background 

0.019 

0.011 

0.007 

0.009 

B° —> K *°7 contamination 

0.003 

0.004 

0.003 

0.002 

Fit bias 

0.008 

- 

- 

0.010 

Total systematic uncertainty 

0.03 

0.05 

0.05 

0.05 

Statistical uncertainty 

0.06 

0.23 

0.22 

0.18 


are obtained from the uncertainties on the parameters describing the shapes and by 
allowing for potential cos 20 and sin 20 modulations. 

To estimate the uncertainties due to the modelling of the PR background the Tl.pr 
parameter is varied between 0.17 and 0.5. The systematic uncertainties related to the 
and A observables depend on the values of the observables themselves: their sizes are 
assessed by varying the damping factor up to 0.5, i.e. reducing the distortions of the 0 
distribution of the PR background by a factor of two compared to the signal ones. For the 
parameter, the systematic uncertainty is estimated by varying R^pr u P t° the fitted 
value obtained for B°—y K*°e + e~. 

The systematic uncertainties from the B° —y K*°j e + e - background are due to the 
uncertainty on the size of the contamination. 

Finally, to estimate possible biases due to the fitting procedure, a large number of 
pseudo-experiments are generated with the number of events observed in data and are 
fitted with the default PDFs. While the A^ 1 and AJjP estimates are not biased, the Ff 
and observables exhibit small biases (less than 10 % of the statistical uncertainties) 
due to the limited size of the data sample and are corrected accordingly. The values of the 
corrections are assigned as uncertainties (labelled as “Fit bias” in Table [2j) . 

The systematic uncertainties are summarised in Table |2j The systematic uncertainties 
on the i%, A^\ Al^ 1 and angular observables in Table l] are obtained by adding 
these contributions in quadrature. They are, in all cases, smaller than the statistical 
uncertainties. 


7.3 Effective q 2 range of the selected B { 


I<*°e + e- 


signal events 


The distribution of the reconstructed q 2 for the signal is obtained using the sPlot tech¬ 
nique |40| based on the B° invariant mass spectrum and shown in Fig. [5j Taking into 
account the effect of event migration in and out the q 2 bin, the average value of the true 
q 2 of the selected signal events is equal to q 2 = 0.17 ± 0.04GeV 2 /c 4 . The acceptance 
as a function of the true q 2 , obtained from the LHCb simulation, is uniform in a large 
domain except close to the limits of the reconstructed q 2 , 0.0004 and 1 GeV 2 /c 4 . Due to 
reconstruction effects, the q 2 effective limits are slightly different. Because of reduced 
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Figure 5: Distribution of the reconstructed q 2 from an sPlot of data (black points). The dashed 
line represents the B° —> K*°e + e~ contribution and the grey area corresponds to the 3.8% 
BK*° 7 e + e - contamination. The solid line is the sum of the two. 


acceptance in the low-g 2 region, the value of the lower q 2 effective limit is increased; because 
of bremsstrahlung radiation, events with a true q 2 greater than 1 GeV 2 /c 4 are accepted by 
the selection and the higher q 2 effective limit is also increased. The values of these effective 
boundaries are obtained by requiring that in the low- and high- q 2 regions the same number 
of events are obtained in a uniform acceptance model and in the LHCb simulation. The 
true q 2 effective region is thus determined to be between 0.002 and 1.12 GeV 2 /c 4 . It is 
checked, using the LHCb simulation, that the average values of the true q 2 and of the 
angular observables evaluated with a uniform acceptance in the region between 0.002 and 
1.12 GeV 2 /c 4 are in agreement with those obtained from the angular fit performed on the 
events selected in the reconstructed q 2 interval 0.0004 to lGeV 2 /c 4 . An uncertainty on 
the q 2 effective limits is assigned as half of the q 2 limit modification. The true q 2 effective 
range is thus from 0.0020 ± 0.0008 to 1.120 ± 0.060 GeV 2 /c 4 . This range should be used to 
compare the A^\ AljP and measurements with predictions. 


8 Summary 

An angular analysis of the B° —>■ K*°e + e~ decay is performed using proton-proton colli¬ 
sion data, corresponding to an integrated luminosity of 3.0 ftT 1 , collected by the LHCb 
experiment in 2011 and 2012. Angular observables are measured for the first time in an 
effective q 2 range from 0.0020 ± 0.0008 to 1.120 ± 0.060 GeV 2 /c 4 . The results are 

F l = 0.16 ±0.06 ±0.03 

4 2) = -0.23 ±0.23 ±0.05 
4 m = ±0.14 ±0.22 ±0.05 
A* e = ±0.10 ±0.18 ±0.05, 
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where the first contribution to the uncertainty is statistical and the second systematic. 
The results are consistent with SM predictions |2,41 j . For the low average value of q 2 of 
this analysis, the formulae relating and A ] f n and C 7 and C' 7 in Eq. [ 3 ] are accurate at 
the 5% level, for SM values of the ratios of Wilson coefficients C 9 /C 7 and C 10 /C 7 . At this 
level of precision and for SM values of C 7i the ratio C 7 /Cy is compatible with zero. This 
determination is more precise than that obtained from the average of the time-dependent 


measurements of CP asymmetry in B° —y K*°(—> decays 16,17 
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